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Abstract: We use a subtraction method to construct NLO corrections in a 
Monte-Carlo event generator for the case of vector boson production in Drell-Yan pro- 
cesses. Our calculations are carried out both for the Bengtsson-Sjostrand-van Zijl (BSZ) 
algorithm and for a modified algorithm proposed by Collins. In the case of the modified 
algorithm, we compute the relation between the parton distribution functions and the 
ones in the MS scheme; this relation is the same as the corresponding relation for DIS. 
For the BSZ algorithm, we show that there is no simple relation. 
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1. Introduction 

In an earlier paper by Collins Jl|, a subtraction method was introduced to consis- 
tently take into account next-to-leading order (NLO) terms for deep-inelastic scatter- 
ing in a Monte-Carlo event generator. In that paper, the method was applied to the 
photon-gluon-fusion process in an event generator that uses the algorithm constructed 
by Bengtsson and Sjostrand (BS) for initial-state showering; such event generators 
are PYTHIA [[§, LEPTO |§ or RAPGAP § x . In this paper, we apply this method to 
massive vector boson production (W or Z production) in hadron-hadron collisions. The 
method contrasts with previous methods, e.g., ||, at incorporating NLO corrections 
by a reweighting of the events generated by showering from the LO matrix elements. 

1 The same methods can be applied to event generators like HERWIG that use a different algorithm, 
but the details of the calculation will be different. 
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The subtraction method is intended to be applicable to all the perturbative parts of 
an event generator, to calculate non-leading corrections to the hard scattering and the 
showering. 

As described in IJ, an event generator using this subtraction method to correct 
the hard scattering generates two classes of events. One class is obtained from the 
LO parton-model process by showering the initial and final state quarks, exactly as 
before. The second class of events is generated by starting with an NLO subprocess and 
showering the partons, again exactly as before, but with one exception, the exception 
being that the hard cross section for the subprocess is equipped with a subtraction 
that correctly compensates the double counting between the two classes of events. 
This removes the part of the NLO term that is included in the combination of the LO 
parton model and showering, and it cancels the singularity in the NLO contribution to 
the cross section. 

We will use two different algorithms for the parton showers. One algorithm is due 
to Bengtsson, Sjdstrand and van Zijl || and is used in PYTHIA and RAPGAP. As 
explained in Ref. the parton density functions (pdf's) to be used in the event gen- 
erator are not those of a standard scheme, but are specific to the showering algorithm. 
We will actually find that there is no simple relation whatever between the parton 
densities to be used with this algorithm and the MS pdf's, and that they are therfore 
also different from the pdf's needed for the corresponding algorithm in DIS. This 
happens because of the changes made by the showering algorithm to the parton model 
kinematics: the changes are different from those in DIS, and correlate the kinematics 
of the partons from each of the incident hadrons. The second algorithm is the one 
defined in [[1J with the specific aim that the relevant parton kinematic variables do not 
get changed by the showering. For this algorithm, we will show that the pdf's for the 
hadron-hadron induced process are the same as in DIS, although both differ from the 
MS definition. 

The NLO corrections for vector boson production come from two different sub- 
processes: gq and qq. In this paper we derive the results for gq subprocesses; this is 
an important correction which is not necessarily suppressed by the factor of a s if the 
gluon density is larger than the quark densities. The generalization to qq subprocess 
will encounter some complications because of the need to treat soft gluon effects. So 
we defer this to the future. 

In Sec. |2], we describe the treatment of vector boson production in an event gener- 
ator. Then we compute the effect of combining the LO cross section with the order a s 
part of the shower; this will be needed as the subtraction term in the NLO calculation. 
We present the calculation for Z bosons, and later, in Sec. |], specify the changes needed 
to treat the production of W bosons. 
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In Sec. H], we carry out the subtraction from the NLO matrix elements to get the 
resulting NLO hard differential cross-sections. The pdf 's in these formulae are not in 
MS scheme, so we show how to relate them to the ordinary MS pdf's. For the new 
algorithm, we find that the relationship between the pdf's is the same as was found with 
deep- inelastic scattering in PJ; while for the BSZ algorithm we show that no simple 
relation appears to be possible. 

Finally, we summarize and discuss our results in Sec. |5|. 



2. Monte- Carlo algorithms 

2.1 QCD improved parton model and vector boson production 

Our calculations are based on the QCD improved parton model, which shows that a 
hard scattering process initiated by two hadrons is the result of an interaction between 
one parton (quark or gluon) from each of the incoming hadrons f8j. 

We consider Z boson production in a collision of two hadrons A and B at a center- 
of-mass energy ^/s. The leading-order subprocess is qq — > Z. As in Fig. [I], the resulting 
cross section da^/dyo for producing a Z boson of rapidity y is obtained by weighting 
the subprocess cross section with the parton distribution functions fi(x a ) and fi{xb), 
and summing over all quark- antiquark combinations in the beam and target: 

r1a LO 

-7^ =Kj2(A* + Vnf i/A (x a ,Ml)fy B (x b ,M* z ). (2.1) 

"I/O i 

The parton momentum fractions are written as 

My My 

Xa = ^\ Xb = -^e- y °. (2.2) 



The parton-level cross-section results in the factor K(Af + V 2 ), where 

K = ^£l T = M = w (2.3) 



and 

Ai = Tf, Vi = Tf -2Q iS m 2 e w (2.4) 

come from the axial and vector couplings in the electro- weak interaction ||. Finally, 
fa/A(£,fJ> 2 ) is the number density of quarks of flavor a in hadron A at fractional mo- 
mentum £ and a renormalization/factorization scale \x. 

In the usual "matrix element" approach, the variable yo is exactly the rapidity of 
the Z boson. This follows from the approximation of giving the incoming partons zero 
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transverse momentum and virtuality. However, in a Monte-Carlo event generator, the 
quark and antiquark are given their correct kinematics. In that case, the variable yo is 
not exactly the rapidity of the Z boson; its precise definition is \ ln(x a /xf>), where x a and 
Xb are the fractional longitudinal momenta of the incoming partons. This is a variable 
that is well-defined in the generation of a particular event, but that is not necessarily 
measurable from the final state of the event. In the context of an event generator, 
Eq. (|2.1| ) must therefore be reinterpreted, not as the lowest-order approximation to the 
physical cross section for a particular Z-boson rapidity, but as the cross section for 
events that inside the program have a certain value for the variable yo- Of course, the 
motivation for using the variable yo is that it approaches the true rapidity y in the limit 
that the parton transverse momenta and virtualities approach zero. 

Since the Z's decay width is small compared to its mass Mz, it is sufficient to 
compute the production cross sections of the bosons. The actually measured cross 
sections of leptons are computed by multiplying the cross sections by the appropriate 
branching ratios. 




Figure 1: Vector boson production in the Drell-Yan model 



With full perturbative QCD corrections taken into account, Eq. ( |2 . 1| ) gets replaced 
by the factorization formula 

^ = £/ fi/Afotffj/B&y)^, (2-5) 

where y is the rapidity of the Z boson, the £'s are the momentum fractions of the 
incoming partons, and da^jdy is a suitably constructed hard scattering cross section. 
Now the sum is over all pairs of parton flavors (quarks and gluons). The formal domain 
of validity of Eq. ( |2.5|) is for the inclusive cross section in the asymptotic 'scaling' limit, 
analogous to the Bjorken limit in DIS, s — > oo with r and y fixed. As is well-known, 
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one effect of higher-order perturbative corrections is the production of vector bosons at 
large transverse momentum (qr)- The discussion above, about the distinction between 
the variables y and yo, alerts us that care will be needed in our application of the 
factorization formula in an event generator. In effect, the quantitative interpretation 
of Eq. ( |2.5| ) gets modified in current MC event generators; notably the lowest-order 
formula Eq. ( |2.1| ) is no longer for do jdy but for da/dy^. Effectively, when only LO hard 
scattering is concerned and when q? is small, the foundation of an event generator's 
algorithm is an appropriate modification of Eq. ( J2.5[ ) that embodies the same physics. 
However, when large gr is concerned, we have to be careful about the difference between 
y and y for all orders of hard scattering, including LO. 

2.2 Parton-shower algorithm 

Now we describe the initial-state shower algorithm for vector boson production used in 
PYTHIA, LEPTO or RAPGAP, as described in §. 

From Pa 3' a' From Pb 




0(Z) 

Figure 2: Initial-state parton shower for vector boson production 

In Fig. |2| is symbolized an example of initial-state parton showering for Z produc- 
tion. The showering algorithm generates partons with certain flavors and virtualities; 
it also generates the splitting variables z 2 n+i, z 2n for each branching, and an azimuthal 
angle for the transverse momentum of each branching. The momentum fraction of each 
space-like line is computed as 

6n+l = T-rn-1* ' 6n = , (2.6) 

llj=0 Z 2i+1 lli=l Z2i 

so that £ 3 = £4 = £2/^2, etc. Since the z's are generated numerically, by an 

algorithm explained below, Eq. (|2.6|) gives the values of the momentum fraction vari- 
ables for all the lines, except for £1 and £2- These values are defined to be x a and 
Xb, which are generated according to the probability distribution corresponding to the 
lowest order cross section, Eq. (|2.1|). 

As in , the "s approach" is used to relate the splitting variables Zi to the parton 
4-momenta. This is done by requiring that 

Sij = &€jS (2.7) 
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both at the hard scattering and at any lower scale in backward showering, where £j and 
£j are of the two resolved partons. This means that the total % has to be increased 
by a factor of \jz in the backward evolution; this defines the relation between z and 
parton kinematics. For instance, in Fig. 0, if line 1 has the highest virtuality, then 

Z\ = {Pl+P2) 2 /{P3+P2) 2 - 

The part of the algorithm used in an event generator that concerns us is as follows: 

1. Generate values of yo and M| from the LO cross section for Z boson production 
Eq. Q2.1p . From these variables, calculate x a and Xb by Eq. 



2. Generate a virtuality Q\ for the incoming quark a, a longitudinal splitting variable 
fraction z\ for the first branching, and an azimuthal angle <j) for this branching. 
The distributions arise from the Sudakov form factor 

*.(*., «?) = -P { - /f " f^ 2 ^ (2-8) 



x 



V- f 1 dZi ^fk(Xa/Zl,Q' 2 ' 
2^ / ~^k-nj\Zx t 



Here, Q^ax * s normally set equal to M§. The Sudakov form factor is the proba- 
bility that the virtuality of quark a is less than Q\. 

3. Iterate the branching for all initial-state and final-state 2 partons until no further 
branchings are possible. 

4. Compute 4- vectors for the momenta of the generated partons. 

2.3 First initial-state branching with BSZ Algorithm 

Later, in Sec. |3|, we will calculate 
the hard-scattering cross section for the 
gg-induced process, Fig. The gluon, 
of momentum p%, comes from hadron A 
and the (anti)-quark, of momentum p%, 
comes from hadron B. The NLO con- 
tribution we want to calculate is to be Figure 3: NLO gluon-quark collision 
accurate when the incoming gluon and 

quark have virtualities and transverse momenta that are small compared with Mz, and 
the intermediate quark, of momentum pi, has a virtuality of order M§ . To avoid double 
counting, it is necessary to subtract the corresponding contribution obtained from the 
showering algorithm applied to the LO partonic cross section, and it is this subtraction 
term that we calculate in this section. 





2 The final-state showering is organized similarly to the initial-state showering. However, we will 
not need it explicitly in this paper. 
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The subtraction is obtained by multiplying the lowest order cross section, from Eq. 
(|2.1|) , by the appropriate part of the showering approximation at order a s . Only the 
gluon-to-quark splitting is relevant for our calculation, and to match with the definition 
of the NLO hard-scattering, it should be calculated when the initial-state partons are 
on-shell and have zero transverse momentum, and the final-state quark is on-shell. 
The showering factor is just the first order term in the expansion of the Sudakov form 
factor in powers of a s (M§). 

The first-order cross section in the showering approximation is 

<i<7 ^ ,[ (Al + Vt)^^C(Ql)P( Zl )U(^Ml)Mx b ,Ml). (2.9) 



dyodQld&dcj) t %, A^Q\ ^ u viy £a 

Here, £3 is the longitudinal momentum fraction of p^, and the splitting kernel is for 
g — > quark + antiquark: P(zi) = P g ^ff(zi) — |(1 — 2zi + 2zf). Because of the way 
in which an event generator uses the lowest-order cross section, it is the variable yo 
that appears in Eq. (|2.9| ) rather than the true rapidity of the Z boson. We will relate 
yo to the true rapidity later. Note that because we are doing a strict expansion in 
powers of a s (M|), the scale argument of the pdf's is M§. The function C{Q\) is a 
cut-off function |l| that gives the maximum value of Qf, and the standard choice is 
C(Qf) = Q{M\ — Ql). We will not discuss other choices of cut-off functions in this 
paper. 

Now we reconstruct the 4- vectors for the momenta q, p±, ps and p[ of the vec- 
tor boson, intermediate quark, the incoming gluon and the outgoing quark. In the 
Bengtsson-Sjostrand-van Zijl's definition ||, they obey the following requirements: 

1. Hadron A is to be moving in the z direction. 

2. Hadron B is to be moving in the —z direction. 

3. The incoming partons, P2 and p% have momentum fractions £2 and £3 relative to 
their parent hadrons, in the sense of light-front components. 

4. p\ = -Ql 

5. p\ = p\ = pf = 0, and P2 and p^ have zero transverse momentum. 

6. q 2 = Ml 6 = x b . 
O1+P2) 2 x a x b s x a 



7. z x 



(P3+P2) 2 £3' 
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In the center-of-mass frame of hadrons A and B, with the components written in 
the order (p°, Pt,P z ), we then have 



Pi 



P2 



u 

Pa 
Pi 



1,0t,-1), 



Pt, 



2v^ 



M 2 +Q\ Q\ 
6 6 



^(1,0, 



Pl-Pl 



1), 



-Pt,P 3 -Pi 



P°2 



Pi,Pt,P Z 2+pI 



with 



Pt 



Qi i 



M 



(2.10) 

(2.11) 

(2.12) 

(2.13) 
(2.14) 

(2.15) 



where is a unit transverse vector in the direction defined by the azimuthal angle <fi. 

Note that after showering, the rapidity of the Z boson is not yo = \ hi(x a /xb) but 
is given by 

yB i h ,4±i! = i ln (Qi±M%, (2.16) 

which does approach m the limit Qi — > 0. 

The natural variables for the LO differential cross section plus the first-order show- 
ering are Qf, yo, and £3. However, they are not so convenient for the NLO corrections. 
So we now transform the cross section in Eq. ( |2.9|) in terms of more convenient variables 
for a hard gluon-quark scattering: y, £ 2 , £3. 

From the above equations, we have 



Ql 



e 2y s& 



M 2 



1 Mz w 
In — - In f 2 , 



which gives the the Jacobian 

d(y ,Ql) 



2e 2 y^ 3 s + M 2 z ] 
(6 + e 2 ^ 2 ) 2 



Then the cross section is 



da. 



(BSZ) 
shower 1 



dy d& d£ 3 ' 



(A 2 + V 2 )^^-C(Ql)P( Zl ) 

J x/ 3 fe,Mi)A fe ,Mi) 2 2f;4 2 y 



(2.17) 
(2.18) 

(2.19) 

(2.20) 
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2.4 First initial-state branching with New algorithm 

In the standard algorithms used for treating parton kinematics, the effect of shower- 
ing changes the relationship between observable quantities and the parton momentum 
fractions from their parton model values. In Ref. [|lj an algorithm was proposed that 
does not suffer from this effect. Here we extend this algorithm to the Drell-Yan process 
by requiring the rapidity y of the Z boson to be the same after we shower the incoming 
partons. (Thus it is not necessary to use a separate variable y to denote | lia(x a /xb).) 

The parton momenta obeys the same requirement as the BSZ algorithm, except 
that items || and [7] are replaced by 

|f q 2 =MJ andy= - In—. 

2 Xb 

Thus the condition on the momentum of p 2 is dropped, and instead the rapidity of 
the Z boson is required to obey the simple parton-model relation to x a and Xb- In 
addition, the first splitting variable z\ is defined by without using the "s" condition. 
Accordingly, the fraction momentum of is required to be £3 = x a jz\ in the new 
algorithm. Note that Eqs. ( ^.10[ )-( p.l(j| ) remain true in the new algorithm. What 
has changed is the relation between the parton kinematics and the variable Xb that 
is generated by the algorithm. In the old algorithm Xb = £2; in the new algorithm 
x b = x a e~ 2y = i z z x e- 2y . 

The first-order cross section in the showering approximation in the new algorithm is 
obtained by transforming the cross section in Eq. ( |2.9| ) to be differential in the variables 
Vi €2, ^3) 4>i by using Eq. ( |2.17| ), which gives the relation between £ 2 and Q\\ 

j (New) / j.j2 

acr shower 1 K \^ (A 2 I V 2 ) ^ 

dydz 2 d£ 3 d<t>~ Y ^ 2 Qi 

^ 2y [i2i,s{2 + j 2 /i,e 2 y) + M 2 ] 

x 77 T 2vc \2 M&> M z)Mxb, M z ), 

where y is the exact rapidity of Z boson. Notice that in this algorithm, although the 
parton density for the parton of momentum p 2 is calculated with momentum fraction 
Xb = e v Mz/\fs, the actual fractional momentum is different (except at zero transverse 
momentum), and is given by a rather complicated formula. 

This difference between the actual parton momentum fraction and the value used 
in the evaluation of the parton density is a characteristic of this algorithm. The utility 
of this apparent inconsistency will become apparent when we compute the relation with 
the results of the "matrix element" method of calculation. 



^ 2 + V 2 )^^C l {Q 2 l )P{z l ) (2.21) 
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3. NLO hard cross section 



3.1 Unsubtracted NLO term 

From standard references (e.g., 0), we find that the unsubtracted cross section asso- 
ciated with gluon-quark scattering corresponding to Fig. |3| is: 

^subtracted i = a s (M 2 ) V2G F M 2 (A 2 + V 2 ) 
d&d&dtdcf) Y 452 3 ( 27r ) ' 

x P + i2+ f 2{iM2z ms, Ml)A(C 2 , Ml). 
st 

The sum is over quark and antiquark flavors, and 



s = (P3+P 2 ) 2 = 66s, (3.2) 
i = -Ql 

u = M 2 Z -t-s = M 2 z + Q\-& £ 3 s. 



From Eq. ( 2.171 ), we have 



dt_ -2e 2 ^2 6[6e 3 g + M|] 
dy (£3 + e 2 ^ 2 ) 2 



(3.3) 



So now we can also write this unsubtracted cross section in terms of the same variables 
as Eq. flT2lD : 



^unsubtracted 1 ^ 2 2 a s (M z ) S 2 + t 2 + 2u M z 

e 2 y(^ 3 s + M 2 



3.2 NLO term with subtraction 



M|)A(6, M|) 



We now subtract the showering term, Eq. ( |2.20| ) or ( |2.21|) , from the above 0(a s ) term. 
A change in the labeling of the parton momentum fractions is in order. The variables 
we have previously used, £ 2 , £3, etc, were tied to a particular structure for the LO 
hard scattering and showering. But for the subtracted NLO cross section we prefer 
something that has the same names for the external partons independently of the 
subprocess considered. So we make the following change of notation: Instead of £3, we 
will write £ a to indicate that a parton from hadron A, and instead of £ 2 , we will write 
£5 to indicate that a parton from hadron B. 
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For the BSZ algorithm: 



1 ' -4nH 
e 2 y[UbS + M 2 z \ \ s 2 + i 2 + 2uM 2 z 



(3.5) 



2d(-t)P^i) 



For the new algorithm: 



o2.y 



-47T 2 t (^a + e 2 ^ 



(s + M: 



(3.6) 



- /,(x 6 , Mf) d(-t) P(^) S f 2 + I e 2 ^ + ^f- 



In both equations s, t, u, and z\ are all functions of £ a and given by the formulae 
earlier, which in terms of the new notation are: 



-u = UbS -M 2 z -Qi, 



3.3 Results for gq subprocess 



(3.7) 

(3.8) 
(3.9) 





Pi P2 



Figure 4: NLO gluon-quark collision 



Now we present the results for the other gluon-quark scattering subprocess, in 
which gluon comes out of hadron B instead of A, from diagrams shown in Fig. [|. The 
formulae for the subtracted cross section are obtained from Eqs. ( |3.5| ) and ( |3.6| ) by 
making the following change: 



2, a 



b, y 



-y- 



(3.10) 
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The resulting subtracted NLO cross sections are: 



For the BSZ algorithm: 



^E + ^^m a ,Ml)f g {t h ,Ml) (3.11) 

-2C 2 (-u)P(z 2 )l. 



dy d£, a d£ b d(f) i -47r 2 u 

e- 2y [^ b s + Mf] / s 2 + u 2 + 2tM| 



fe + e- 2 ^a) 2 * 2 



For the new algorithm: 

, (New) / A ^-2 \ _-2y 

hard2 " ^E (A 2 + O ^¥ 7^-^7-2 Mf) 



dydZ a d£ b d<P i l ' -4tt 2 m (£ 6 + e -%£ a 

s 2 + M 2 + 2tM| 
s 



x\m,Mt) ° T ";" J " z (i + Ml) (3.12) 



/A M 2 
- M|) C 2 (-«) P(z 2 ) S 2 + £ e- 2 y + 



Here, z 2 = x b /£ b . 



3.4 Comparison with MS scheme 

As discussed in ffl, after we have obtained the NLO corrections, it is necessary to 
find the relation between the scheme of the pdf's used in the event generator and the 
commonly used MS scheme. 3 As in [|IJ we will compare the same cross section computed 
with the MS pdf's and with the formulae used in the Monte-Carlo approach. We will 
choose to calculate da/dy, since 

• This cross section can be calculated analytically both in the normal factorization 
method and from the algorithm used in the event generator. (The inclusive cross 
section allows an integral over the showering probabilities.) 

3 The need for the change of scheme is not apparent in other work on merging parton showers and 
matrix elements || . Our work relies on a deeper analysis of the derivation of the whole algorithm used 
in a Monte-Carlo event generator, rather than just a consideration of the normalization of a particular 
inclusive cross section. It is only in this context that it becomes apparent that the MS parton densities 
are not the appropriate ones. The parton densities used in an event generator implicitly include 
observed jets that result from initial-state showering, and therefore the definition of the parton density 
is directly tied to the definition of the showering algorithm. One way to ensure that the pdf's for the 
event generator are the same as MS pdf's is to adjust the cut off function C suitably as is done by 
Potter §. 
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• The lowest order has a factorized dependence on both parton momentum frac- 
tions, so there is exactly enough information to extract the parton densities. 

As already discussed in |]], the unsubtracted part of the cross section is the same 
in both methods of calculation, so that the result of the calculation depends only on 
the subtraction terms. The relation between the pdf's is therefore independent of the 
details of the electro-weak couplings for example. 

It is less obvious that the relation between the pdf's could be unaffected by the 
differences between the kinematics for the DIS and DY processes. This is what we will 
find, nevertheless. 

In order to get a complete determination of the parton densities, we consider a 
modified cross section, in which only a single quark and antiquark flavor annihilate, 
and in which Mz is given an arbitrary value. Moreover, the electroweak couplings form 
a common factor, and it is convenient to remove them and to arrange that the lowest 
order cross section is just fi(x a )fi(xb); this we call the structure function !Fi(y, M|). 



First we take the formula for the structure function in the MS scheme ||10||: 



Hy,M z ) = ri MS \x a ,^)fr s \x b ,^) 

Z7T JXa JXh 



X 



ff% b , 5(6 - x b ) i (p{ Zi ) In 2fet ^ .. Xb) +z x {l- Zl ) 

Sa \ Xb\<,a T X a j 



with 



and 



Kb- %b J ?6 - x b 

+ (i <-> i, a <-» b, y <-> —y) > + first-order quark terms + 0{a 2 s ) (3.13) 



r c f t <■ \ _ x b( T + 66) V 2 + ( r ~ 66) 2 ] (o-, A \ 
C(U6J- Q^^ Xb + ^ Xa )^ + Xb ) ' ^ 

rrcfc m T Jj + 6 6) [6 6 2 x a + r(£ a x b + 2 6 x a )] 



Here the variables x a and x b are defined to be exactly the "parton-model" values x a = 
\fre y and x b = y/re~ y , with r = M|/s. The variable Zi in the splitting kernel is 
defined as Z\ = x a /£ a , the same as in the new algorithm. 

This structure function must equal the same structure function given by the Monte- 
Carlo calculation. Since the expressions from the new algorithm are simpler and have 
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similar structure, we will first compare with the structure obtained using the new 
algorithm . 

We start from Eq. fl3.6| ) for the cross section. Using the functions defined in Eq. 
( 3.14Q and Eq. ( 3.15|) , we rewrite so as to obtain a formula with a similar structure to 
the MS formula: 



2tt 



1 di a f di b {f^\i ai M 2 z )fi 



(x b , Ml 



x5(^ b -x b )C 1 (Q 2 1 )P(z 1 



2l r„, ^ e b (e a -xi)-(T+^ b ) 



X a £a (£ft % b ) (£a,X b + £,bX a ) 



+ft w) (^M^fr w> (^M 2 



2 \ r (new) . 



+ {i <-> i, a <-> b, y <-> — y) | 
+first-order quark terms + O(a^) 



+ # C (£a,& 



(3.16) 



Comparison of Eq. (3.16) and Eq. (3.13) shows that 



x/i new )(x,Mf) 



x/i MS) (^,Mi) 



a s (M. 



2tt 



rfe|4 MS) ce, Af) w - «) + z(i - z)] 



+first-order quark terms + O(af), 



(3.17) 



and 2 = the same as in Ref. |T[. 

When we attempt a similar calculation with the BSZ algorithm, we find that the 
parton densities in the subtraction term in the equivalent of Eq. (|3.16|) are no longer a 
product of fg(£ a ) an d fj{x b ). Instead the argument, x b = y/re~ y , of the second parton 
density is replaced by a complicated function of the kinematic variables. This does 
not match the structure of the MS formula Eq. (|3.13|) , and it is not possible in any 
simple way to extract a relation between the parton densities for the BSZ algorithm 
and for the MS scheme. The problem is that the showering of the parton on the A 
side has affected the kinematics of the parton on the B side and that it is the modified 
kinematics that are used in the corresponding parton density. As far as we can see, a 
correct analysis can only be done by investigating the problem in terms of unintegrated 
parton densities. Since this is a much more complicated problem, we shall not attempt 
it here. 
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4. W production 



In the above sections we described only Z production since the formulae are for anni- 
hilation of a quark % with its antiquark %. For W production the results are similar, 
only that one needs a different flavor of antiquark ll and an appropriate change in the 
overall coupling: 

V2G F (A 2 + V 2 ) M\ — V2G F \V lV \ 2 M 2 w , (4.1) 
where Vw is an element of CKM matrix. 



5. Conclusion 

We showed how to incorporate the gluon-quark processes in a Monte-Carlo event gen- 
erator using the subtraction method proposed by Collins in Ref. |IJ. We also analyzed 
the exact parton kinematics used in the BSZ algorithm, and observed that the factor- 
ization theorem for the Drell-Yan process is used in a modified form compared with 
the form normally used for the inclusive cross section da/dy. 

When we computed the relation between the parton densities for the event gener- 
ator and the MS densities, we found the same relation as found in |I| in the context of 
deep-inelastic scattering, but only if we used the new algorithm for parton kinematics 
that was proposed in This supports the hypothesis of process independence of 
the pdf 's. For the case of the normal BSZ algorithm, we found that the effect of the 
correlated parton kinematics appears not to permit us to obtain a simple relation. 
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